Antispiral waves are sources in oscillatory reaction-diffusion media * 
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Spiral and antispiral waves are studied numerically in two examples of oscillatory reaction-diffusion 
media and analytically in the corresponding complex Ginzburg-Landau equation (CGLE). We argue 
that both these structures are sources of waves in oscillatory media, which are distinguished only by 
the sign of the phase velocity of the emitted waves. Using known analytical results in the CGLE, we 
obtain a criterion for the CGLE coefficients that predicts whether antispirals or spirals will occur in 
the corresponding reaction-diffusion systems. We apply this criterion to the FitzHugh-Nagumo and 
Brusselator models by deriving the CGLE near the Hopf bifurcations of the respective equations. 
Numerical simulations of the full reaction-diffusion equations confirm the validity of our simple 
criterion near the onset of oscillations. They also reveal that antispirals often occur near the onset 
and turn into spirals further away from it. The transition from antispirals to spirals is characterized 
by a divergence in the wavelength. A tentative interpretaion of recent experimental observations of 
antispiral waves in the Belousov-Zhabotinsky reaction in a microemulsion is given. 
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I. INTRODUCTION 

Chemical pattern formation results typically from the 
interplay of reaction and diffusion and occurs in many 
reactions in solutions and gels Q and in heterogeneous 
catalysis 0- Reaction-diffusion processes are also be- 
lieved to be at the heart of morphogenesis in biological 
systems 0, Rotating spiral waves are probably the 
most typical structures investigated so far. 

They have been found initially in the Belousov- 
Zhabotinsky reaction Since then, they have been 
frequently observed in a variety of experimental setups 
including aggregation of slime molds [3, catalytic CO 
oxidation 0| , cardiac tissue Q , and intracellular calcium 
dynamics in frog eggs 0, and glycolytic activity in ex- 
tracts of yeast cells [l^l ■ 

Spirals that organize the surrounding medium by reg- 
ular emission of waves are well established under both 
oscillatory and excitable conditions. If several spirals 
populate the system then the waves emitted by neigh- 
bouring spirals annihilate each other at the boundary 
of the spirals's spatial domains. Such boundaries con- 
stitute wave sinks. Intuition tell us that these passive 
objects (where waves arrive and die) do not influence the 
medium beyond its closest surrounding. 

Recently, Vanag and Epstein reported on experiments 
of the oscillatory Belousov-Zhabotinsky reaction in an 
Aerosol OT microemulsion (BZ-AOT) In their quasi 
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two-dimensional system chemical concentration patterns 
are arranged as so-called "antispirals" or "antitargets" . 
Typical snapshots of the system display closely packed 
domains each containing a single antispiral. As time pro- 
gresses the concentration waves propagate towards the 
center of the respective antispiral. Such inward propaga- 
tion seemingly contradicts our intuition of spiral cores as 
wave sources. A word of caution towards the terminol- 
ogy of spirals and antispirals is in place here: Spirals and 
antispirals are not structures that can appear simultane- 
ously for given parameters. Instead, they occur in differ- 
ent complementary regions of the parameter space. The 
terminology of antispirals used here is rather recent and 
occurs only after the cited experimental paper, that is 
more careful and talks of inwardly rotating spiral waves. 
It should also be stressed that many authors have found 
(what is now often called) antispirals or antitargets in 
numerical simulations in continuous and discrete oscil- 
latory and chaotic media and simply classified them as 
spirals or targets see e.g. [l^. ITa. Ill 111 111 ITtI IT8l|. 
More recent reports have adapted the new name antispi- 
ral in the sense it is used here 0, 153 . The distinction 
achieved here by the prefix anti is sometimes also referred 
to as positive dispersion (spiral) and negative dispersion 
(antispiral) [T^ . 

Since, the phenomena of antispirals and antitargets oc- 
cur frequently in the complex Ginzburg-Landau equation 
|2ll | , they should be generic in oscillatory media near on- 
set. What comes as a surprise is rather that over so 
many years only spirals have been observed and the first 
experimental report of an antispiral occured only very 
recently |lll | . This experiment lead to a discussion about 
the "mechanism" of antispirals. Vanag and Epstein sug- 
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gested that the phenomenon may be related to a Hopf 
bifurcation with finite wavenumber (wave bifurcation) in 
a extended Oregonator model, where they find group 
and phase velocities with opposite signs. Simulations 
and experiments show also "packet waves" , which are 
wave groups that move in opposite direction to the mo- 
tion of the constituting waves . While this argument 
seems valid for the quasi one-dimensional packet waves, 
it does, in our opinion, not apply to the two-dimensional 
case under consideration. After all, near a wave bifurca- 
tion in a two-dimensional medium, one typically observes 
traveling or standing stripes or hexagons and not spiral 
waves. Earlier on, Nicola et al. have reported related 
wave groups near a Turing-wave bifurcation and named 
them "drifting pattern domains" [23l |. 

Gong and Christini have recently investigated 
the CGLE and prototypical two-component reaction- 
diffusion (RD) models and conjectured that antispirals 
only occur near the onset of oscillations |23|. We will 
provide further evidence for their claim. They also sug- 
gest an analytical argument and criterion for the appear- 
ance of antispirals in the CGLE coefficients. Recently, 
we have shown that the criterion for the occurrence of 
antispirals in RD models has to be modified and that an- 
tispirals in the corresponding CGLE may turn out to be 
spirals in the original RD model |2^. Here we present 
the complete analytical derivation of the differing criteria 
for antispirals in the CGLE and in the corresponding RD 
models. 

When studying dynamical phenomena in RD mod- 
els near the onset of pattern formation then the anal- 
ysis of the cor resp onding amplitude equation can be in- 
structive in [la I23. The CGLE represents the am- 
plitude equation for pattern dynamics near the onset of 
oscillations via supercritical Hopf bifurcations [IJ. The 
CGLE for arbitrary RD systems is obtained following 
exact transformation rules 0, ll^. Here, we validate 
our criterion for antispirals in RD systems for two simple 
reaction-diffusion (RD) systems: the FitzHugh-Nagumo 
and Brusselator models by deriving the CGLE coeffi- 
cients for these models, applying our criterion for anti- 
spirals and comparing the result to direct numerical sim- 
ulations near a Hopf bifurcation. In the CGLE, domains 
of the parameter space can be classified according to the 
relative directions (signs) of phase- and group-velocities 
in rotating waves (spirals or antispirals) in line with ar- 
guments which Y. Kuramoto and T. Tsuzuki used first to 
analyse wave sources of the Kuramoto-Sivashinsky equa- 
tion 1^3 ■ Boundaries of the parameter domains are given 
by zero phase or group velocities which yields analyti- 
cal criteria. Extensive numerical simulations of both RD 
models analysed here, corroborate this criterion's predic- 
tion near the onset of oscillations and extend it when the 
model's parameters are driven away form threshold. The 
simulations also suggest that antispirals typically disap- 
pear far away from the onset of oscillations. 

This paper is organised as follows: In the next sec- 
tion we will discuss theoretically the distinction between 



spirals and antispirals within the CGLE framework and 
analytically derive the criterion for either occurrence in 
RD models near the onset of oscillations. In the third sec- 
tion we explore spirals and antispirals in two RD models 
near and far from the onset of oscillations. We will end 
this paper with a short summary and discussion of the 
main results. 



II. WEAKLY NONLINEAR THEORY OF 
SPIRALS AND ANTISPIRALS 

A. RD systems near a oscillatory bifurcation 

A general reaction-diffusion system in two dimensions 
may be described by the partial differential equation 



(1) 



where u{x,i) is a vector of space- and time-dependent 
concentrations, / is a nonlinear vectorial function, V is 
a diffusion tensor and /i is a control parameter. Realistic 
models of chemical pattern forming systems have been 
proposed in this form including recent models of the BZ- 
AOT system by Vanag et al. |^ and Yang et al. 

■ These include many reaction species, respectively 
components of u, and complex reaction terms /. How- 
ever, in this paper we shall study simple models namely 
the Brusselator [23| or the FitzHugh-Nagumo model {3^ . 
(Anti)spiral waves in the latter two models will be stud- 
ied in Sec. IIIII and explained by results of the weakly 
nonlinear analysis that we carry out in the present sec- 
tion. 

Near the supercritical onset of homogeneous oscilla- 
tions in the RD system with frequency fl and eigenvector 
Ui the vector of concentrations u may be decomposed as 



u{x, t) ~ Uq + UiA{x, t) e -I- u\A*{x, t) e 



-int 



(2) 



The evolution of modulations A{x,t) — y^A{x,t)e^''"* of 
a homogeneous oscillation is described by the CGLE 



dtA = A + il + ici)AA - (1 - ic3)\A\'^A 



(3) 



Here e = (/i — /ic)//^c measures the distance from the 
threshold /ic, ci (03) give the linear (nonlinear) disper- 
sion and Co an overall linear frequency shift. The coor- 
dinates X = x/e^ and t = e- of the CGLE are defined 
by characteristic spatial and temporal scales ^, r and de- 
scribe slow modulations in space and time. Note, that 
intrinsic CGLE frequencies uj result only in a small cor- 
rection of order euj to the original frequency of the system 
at the Hopf bifurcation threshold. In the following, vari- 
ables with (without) tilde will be used in the original RD 
system (the derived CGLE). 

Numerical simulations of the CGLE in a two- 
dimensional system with zero flux boundary conditions 
provide examples of antispiral waves as shown in Fig.^ 
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FIG. 1; Numerical simulation of the CGLE revealing antispirals and antitargets for ci — 1,C3 = 0.5. In (a) the real part 
of A, in a system of size 250*250 with zero flux boundary conditions and after a transient of t = 125, is shown. White (black) 
areas correspond to maximum (minimum) values. The black arrows indicate the direction of propagation of phase waves and 
the white arrow denotes the rotation of an individual spiral. The inwardly propagating target waves are induced by a small 
oscillating heterogeneity in their center. The domain boundaries, where new waves emerge and split, are visible between the 
antispirals. In (b) a space-time plot of 7?e[yl], along a horizontal cut at y = 60 (this line contains the cores of the antispirals at 
X = 75 and 175), is shown. After a transient of t = 125 the final state leads to the configuration shown in (a). Initially the cores 
are placed in the homogeneously oscillating background. This homogeneous background is suppressed as the antispirals grow 
until the whole system is filled by the antispirals at f = 60. Note how new waves emerge and split at the domain boundary 
at a: = 125. The evolution of a small perturbation in (b) illustrates the different sign of the phase and group velocities Vph 
and Vgr (the artificial perturbation is applied at a: = 90, t = 65 and propagates with I'gr). The space-time plot shown in (c) 
corresponds to the same situation than in (b) but a fast homogeneous oscillation with SI — 0.02 has been added to the phase 
of A. This illustrates the dynamics in the underlying RD system which, for the value of Q, chosen, also shows antispirals (i.e. 
the phase velocity Vph and the group velocity Vgr have the same direction as in the CGLE). 



B. Spirals and antispirals in the CGLE 

These (anti)spiral waves of the two-dimensional CGLE 
have the form |35| 

^(r,0,t) = F(r)c^('^^ + /(^'*)) , (4) 

in polar coordinates (r, 6). The asymptotic behaviour for 
r ^ is df{r,t)/dr ^ r,F{r) ^ r and denotes a topo- 
logical defect 1^1 = of charge a in the spiral core. We 
will focus on one-armed spiral waves with a = +1(— 1) 
which are related by mirror symmetry and rotate clock- 
wise (counterclockwise). Other authors term these so- 
lutions spiral and antispiral [3^ (we will not adopt this 
terminology here). In this paper, both left-handed or 
right-handed structures can be spirals (outward propaga- 
tion of phase waves) or antispirals (inward propagation 
of phase waves) depending on parameters. 

Asymptotically for r — > oo these phase waves are char- 
acterised by F{r) ^ \/T~q^ and f{r,t) g^r - ujs{qs)t 
with a selected wavenumber qs and the corresponding 
frequency ujsiqs) = — C3 + ^^(ci -|- C3). The selected 



wavenumber Qs is a function of the parameters ci and 
C3 given by the nonlinear eigenvalue problem resulting 
form inserting Eq. igj into the CGLE ^ |35j| . 

The asymptotic wave field /(r, t) may be rewritten as 
qs{r — Vpht) with the phase velocity Vph = ^s/Qs ■ The 
propagation of small perturbations of the wave is de- 
scribed by the group velocity Vgr = dujs/dqs = 2qs(ci -t- 
C3) . For positive group velocity the spiral or antispiral 
acts as a source that organises the surrounding pattern. 
Phase and group velocity do not necessarily point into 
the same direction. The phase waves move outward (in- 
ward) for positive (negative) phase velocity Vph in the 
CGLE coordinates or Vph in the RD system, respectively. 
Hence, the signs of the selected group and phase veloci- 
ties constitute the defining quantities for the occurrence 
of spirals or antispirals. 

In the following we calculate the parameter dependence 
of the introduced velocities and focus on their signs. The 
selected wavenumber gs(ci, C3) needs to be computed nu- 
merically in g eneral, only its asymptotic are known an- 
alytically [33 ■ However, for the one-dimensional (ID) 
analog of the (anti)spiral wave (i.e. for the ID CGLE 
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in [0, oo) with Dirichlct boundary condition: A{0, t) 
0) the corresponding selected wavenumbcr qsi{ci,cs) 



known analytically |35ll3 



IS 



where a^ci^cs) 



3a(ci,C3)+ /9a(ci,C3)2 cg + 2ciQ!(ci, cg)^ 



2(ci +C3) - V 4(ci + C3)2 



Cl + C3 



for Cl 



4(-2c3 + 9c3 + 7ci) 



C3 ^ 



/3ci(8(ci - C3)2 + 9(1 + ciC3)2 - 4ciC3)i/2 + ci(5 - 9ciC3) - 4c3 



(5a) 
(5b) 



The function 01(01,03) is symmetric and 9^1(01,03) anti- 
symmetric under the substitution (ci — > — Ci,C3 — > 
— C3). The same property holds for the two-dimensional 
(anti)spirals as can be deduced from the form of the non- 
linear eigenvalue problem [35l |: 

9s(ci,C3) = -g,,(-ci, -C3), (6a) 

qsici^cs) = , forci+C3=0, (6b) 

gs(ci,C3) > , forci+C3>0, (6c) 

9s(ci,C3) < , forci+C3<0. (6d) 

Hence, Qs takes the same sign as ci -I-C3. The correspond- 
ing frequency Ws(ci,C3) in the CGLE fulfils ujsicijCs) = 
—uj.^(—ci, — C3) as had already been noted by PauUet et 

The curve where the selected frequency lOsi vanishes 
has been computed numerically in the (ci, C3) parameter 
space {i.e. we computed the parameters ci and C3 for 
which LOsiilsiici, C3)) = 0). It turns out to be different 
from the diagonal C1+C3 = where qs = 0. Hence, within 
the CGLE the parameter space, we can distinguish four 
different domains, as shown in Fig. |21 When crossing 
the boundaries between these domains, the phase veloc- 
ity Vph in the CGLE switches sign between outward and 
inward. Note, that the spiral wave length and Vph di- 
verge on the diagonal and Vph crosses zero where ujs ~ 0. 
In the coordinate frame of the CGLE, there are two do- 
mains with spirals and two with antispirals. This result 
is in agreement with that of Gong et al. and addi- 
tionally covers two quadrants previously not considered, 
where the transition curve Vph = lies. Inserting Eq. H6a|) 
into Vgr ~ 2qs{c\ +C3) one finds that the group velocities 
Vgr of both spirals and antispirals never become negative. 



C. Spirals and antispirals in RD systems 

Now wc return to the initial RD system. Therein the 
asymptotic concentration waves u ~ Q^g^i^-'^pht) -^^itli 
wavenumbcr qs = \/£qs/(, read, after inserting Eq. Q 
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FIG. 2: Parameter space (01,03) of the complex Ginzburg- 
Landau equation (CGLE) and domains where spiral or anti- 
spiral wave solutions exist in the corresponding RD systems 
(where we assumed that SI > 0). The group velocity of spi- 
rals and antispirals is non-negative everywhere. Note that in 
the CGLE the phase velocity vanishes along the dashed curve 
and antispirals (spirals) occur in the upper (lower) shaded and 
lower (upper) white regions. The dot-dashed curves denote 
the Benjamin- Feir instability (1 -I- C1C3 = 0). The white circle 
indicates the location of parameter values used in Fig. 



u{f,9,i) = ito + Mi\/eWl 



X 



^i[a9 + qsf- - {-Vl + e{u}s - co))i] 
+C.C. , (7) 



which yields 



Vph 



-il + -{uis - Co) /qs 

T 



(8) 
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The description given by the CGLE (O is accurate 
only if e 1. In this limit Vph may be simplified as 
Vph ~ ^fl/qs ^ ^fl/qs (note that the same result fol- 
lows for the complex conjugate (c.c.) in Eq. O, where 
two negative signs cancel). For the group velocity we 

find Vgr = Vgr^^^ ~ Vgr which does never get negative. 
Therefore, both spirals and antispirals will act as organ- 
ising centers of the surrounding concentration pattern. 

In order to determine the parameter dependence of the 
phase velocity Vph ~ —fl/qs, we first assume 51 > 0. Then 
spirals with positive phase velocity occur in RD systems 
with corresponding ci -f C3 < as follows from Eq. (|Sc|l . 
Antispirals have qg > and arise in RD systems for which 
ci + C3 > 0. Consequently, spirals and antispirals can 
never occur simultaneously within a single homogeneous 
system. These analytical results are summarised in Fig.|2 
(where 51 > was chosen). 

Let us now briefly comment on the case of 51 < (note 
that the remarks that follow are not important since by 
convention a 51 > is considered in order to calculate the 
amplitude equation ©). Assuming the opposite sign of 
the primary oscillation 51 ^ — 51, yields the same results 
for the group and phase velocities. For negative 51 we find 
ci — !■ — ci,C3 — C3,a; —uj as follows from complex 
conjugation of Eq. |j2J) and qg —qs as in Eq. (|^ . 
This gives the condition ci + C3 < for the existence of 
antispirals if 51 < 0. 

Using Eqs. (jS)) and (|Hc|l . we compare the temporal fre- 
quencies / = —qsVph measured for a spiral (/g) or anti- 
spiral (/as) wave with that of the homogeneous oscilla- 
tion /bulk 

/bulk = ^ + -{co + cz) , (9a) 

T 

/S = 51 + -(CO + C3) + -ql\ci + C3I , (9b) 
r r 

/as = + -(co + C3) - -ql\ci + C3I . (9c) 

T T 

We conclude /g > /b^j^ / a c; < /^uik in agreement 
with experimental observations [ij • 

Let us finish this section with a short summary. The 
observed antispiral waves with phase velocity pointing 
inward were shown to be organising centers since their 
group velocity points outward. Antispirals in oscillatory 
media therefore determine the surrounding pattern as it 
is the case for spiral waves. We found that a spiral in 
the CGLE coordinate frame may represent an antispiral 
in the original RD system and vice versa (see Fig. [SJl. 
Altogether, in RD systems near the onset of oscillations, 
the antispirals are predicted to occur if the corresponding 
CGLE coefficients fulfil ci -f- C3 > 0. 

III. REACTION-DIFFUSION MODELS 

In order to test the general predictions of the previous 
section we will now address our attention to a specific 
RD systems. The RD systems of activator-inhibitor type 



are defined by 

dill = /(u,w) + V^u , (10a) 
= g{u,v) + SV^v , (10b) 

where S is the ratio of diffusion constants and the func- 
tions f{u,v) and g(u,v) define the dynamics of the acti- 
vator u(x,t) and inhibitor v{x,t), respectively. 

In the following we will consider two prototypical ex- 
amples of activator-inhibitor dynamics widely analysed 
in the literature: the FitzHugh-Nagumo and Brussela- 
tor models. In both cases we will proceed as follows: 
First we will study their fixed points and analyse their 
linear stability. The next step will be the derivation of 
the coefficients ci and C3 of the CGLE They can be 
used to decide the question whether spirals or antispirals 
appear near the Hopf bifurcation. Finally we will com- 
pare these prediction with numerical simulations of both 
models near and far from onset. 



A. FitzHugh-Nagumo model 

The FitzHugh-Nagumo (FHN) dynamics !33| is defined 

f{u,v) = u - — -V , (11a) 
g{u,v) = e{u--iv + (3), (lib) 

where e > is the ratio between the time scales of both 
fields, and (3 and 7 are parameters that determine the 
number of fixed points. The coordinates of these fixed 
points are independent of e and are given by the roots of 
a cubic polynomial. In the following we will only consider 
7 > and restrict our analysis to the case where a unique 
fixed point {uq,vq) exists. This is the case, for example, 
for 7=5 where the coordinates of the fixed point are 
uo - (( + 9/32 _ 3/3)i - 1)7(^1 + 9/32 - 3/3)i and 

vo = (wo +/3)/7- 

Choosing e as the control parameter and keep- 
ing /3, 5 and 7 constant, a linear stability anal- 
ysis of (uo,wo) shows that it is unstable to pe- 
riodic oscillations if £ < e^, where = 
(1 — Uo)/7. This Hopf bifurcation has a frequency 51 = 
■\/(l — Wg)/7 — (1 — Uq)2 at onset. The fixed point may 
also become unstable to spatially periodic perturba- 
tions with wavenumber qT ~ \/ ((1 — Uq)(5 — 7^)/ (2(5) 
(i.e. a Turing instability) if e < e^, where = 

((2 - (1 - ul)^ + 2Vl -(l-«g)7)'5)/7'. 

These two instabilities may occur simultaneously at a 
codimension-2 instability. The thick full line in Fig.|3fa) 
shows the location of this codimcnsion-2 instability in the 
parameter space (/3, yfS) for 7 = | • In the following we 
will restrict our analysis to the Hopf instability (ie the 
region below the codimension-2 line in Fig.|3Ia)). 

Using the methods explained in H^H^, the derivation 
of the coefficients c\ and C3 is straightforward. Here we 
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FIG. 3: The gray areas show the parameter space regions where antispirals are predicted to exist near the Hopf instabihty 
threshold for the FHN and Brusselator models (in (a) and (b), respectively). The thick full lines separates the regions where 
either the Turing or Hopf instability appears first as the control parameter e (parameter 6) is decreased (increased) away from 
threshold for the FHN model in (a) (respectively, in (b) the Brusselator model). Inside the dot-dashed lines the Benjamin- Feir 
instability is predicted. The thin dotted lines indicate the range of parameter values considered in the numerical simulations 
of (anti)spirals summarised in Figs. 01 and |n] The Hopf bifurcation considered in both figures is supercritical (see main text). 
In (a) we consider the parameter space (/3, V^) of the FHN model (Eqs. illl ) with 7 = |. For /3 > 2/3 the fixed point {uo, vq) 
is stable and, if e is small enough, excitability is possible. In (b) the control parameter (a, \/5) of Eq. 11211 is plotted. 



will only quote the main results of this derivation. The 
linear dispersion coefficient is given by 



Note that this coefficient vanishes if both fields diffuse 
virith equal strength. If (1 — 7 — 7Mg)/(l — 7 + juq) > 
the Hopf bifurcation is supercritical, otherwise it is sub- 
critical and the CGLE can not be applied. For 7=5 
(the case shown in Fig. Ofa)) the Hopf instability is su- 
percritical for any value of f3. The result for the nonlinear 
dispersion coefficient C3 is more complicated: 

_ 3~3j-7ul + 3j4 I 7 

We may now use the expression for the coefficients Ci 
and C3 to evaluate the conditions for the spiral-antispiral 
transition ci -I- C3 = and the Benjamin- Feir instability 
1 + C1C3 = in the parameter space (/3, v^). These two 
lines are plotted in Fig.l^ta). Note that the region where 
antispirals are predicted to occur near threshold, exists 
only for small values of (3 (i.e. far from the excitability 
region) and vanishes if S is too big or if 5 ^ 0. More- 
over, for /3 = 0, where the model exhibits the symmetry 
(u, v) {—u, — w), only antispirals are predicted. 

The predicted spiral-antispiral transition at ci -I- C3 = 
is strictly valid only close to supercritical Hopf bifurca- 
tions {i.e. for e ^ ^h)- In order to test this and in addi- 
tion to investigate the behaviour far from threshold, we 



have studied the FHN model numerically. We have done 
extensive numerical simulations of lD-(anti)spiral ana- 
logues in Eq. (|ll|l near and far from threshold. These 
simulations were performed in very long systems (typi- 
cally including hundreds of spiral wavelengths) with the 
following boundary conditions: in the right side a zero 
flux condition is imposed and in the left side the field is 
kept to u = mq and v = vq. The output of these nu- 
merical simulations, see Fig. 0] is insensitive to initial 
conditions. If we wait long enough, the boundary condi- 
tion in the left side (the "core") will select phase waves 
with a particular Qsi which eventually invade the rest of 
the system (since they always have positive Vgr). Antispi- 
rals (spirals) will have negative (positive) phase velocity 
Vph- Near the Hopf threshold the observed behaviour is 
equivalent to the one predicted by the CGLE. But, as 
e is decreased, a transition from antispirals to spirals is 
seen. This transition is related to a change in the sign 
of qsi and is not captured by the CGLE description (see 
below) . 

In order to assess the previous results for the ID- 
(anti)spiral analog we have also carried out a number 
of 2D simulations of the FHN model near onset. The 
outcome of these simulations is shown in Fig. ^b) and 
confirms that, at least near the onset, the ID analog 
provides a good description of the spiral-antispiral tran- 
sition mechanism in the 2D system. To see this more 
quantitatively, in Fig. [S] we plot, against e, the selected 
wavenumbers g^, g^i, in ID and 2D, and the correspond- 
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FIG. 4: Summary of the results of numerical simulations of ID- and 2D-(anti)spirals in the FHN model. The thick full line 
in (a) indicates the location of the Hopf instability threshold. The white circle indicates the spiral-antispiral transition point 
at onset, as predicted by the the criterion ci + cs = (c/. Fig. I^J. The dashed curve indicates the spiral-antispiral transition 
of the lD-(anti)spiral analog which was computed numerically. The thick dotted line showd the antispiral-spiral transition 
occurring for small values of e. In (b) a magnification of the region where the spiral-antispiral transition occurs near onset is 
shown. The squares and diamonds show the location of 2D numerical simulations (c/. Fig.0, resulting in spirals and antispiral 
respectively. The parameters are V5 = 1 and 7 = i. Refer to Fig. 01 for a resume of the location of these simulations in the 
parameter spaces (/3, . 



ing wavcnumbers in the CGLE in ID and 2D (using the 
coefficients derived in this section), for the case /3 = 
and (5 = 1. Note, in the inset, the dramatic disagreement 
between the measured wavenumbers and the predicted by 
the CGLE as the distance to the threshold is increased. 
This figure makes explicit the change of sign of qa occur- 
ring for small values of e. Also note that, although the 
measured qs gets apart from qsi as the distance to onset 
increases, the qualitative behaviour is the same. 

B. Brusselator model 

Let us finally address the Brusselator dynamics [ssf. 
This model is defined by: 

f{u,v) = a - (& + 1)m + u^w , (12a) 
g{u,v) ~ bu — u^v, (12b) 

where a and b are two parameters (assumed in the fol- 
lowing to be positive). This model has a unique fixed 
point (uq, vq) = (a, b/a), for any value of the parameters 
a, b and S. 

Similarly to the FitzHugh-Nagumo model, we may now 
examine the stability of the fixed point. Let us take b as 
the control parameter. Depending on the values of the 
parameters a and (5, the fixed point (uo,vo) may become 
unstable in two different ways. If 6 > 6^, where 6^ = 1 + 
a^, then a Hopf instability with frequency fl = a occurs. 



In the other hand, if 6 > 6^, where b^ = (1 + a/VS)"^, 
a Turing instability with wavenumber q^ = a/v^ takes 
place. These instabilities may occur simultaneously in 
a codimension-2 line. This line and the regions in the 
parameter space (a, VS) where either the Hopf or Turing 
instabilities appear first, as the control parameter 6 is 
increased, are shown in Fig.|3Jb). 

A short calculation shows that the coefficients ci 
and C3 are 

(5-1 , 7a2 - 4 - 4a4 

and ^ ^-^^^ . 

The spiral/ antispiral transition line ci+c-^ = is plotted 
in the parameter space (a, v^) in Fig.OJb). Note that the 
region where antispirals are predicted is rather small and 
is located near the area where of the Turing instability 
appears first. Also in this case the antispirals will occur 
only for intermediate values of the diffusion ratio S. 

We have performed numerical simulations also for this 
model. As for the FHN model, we investigated ID-spirals 
and antispirals (with the same boundary conditions as in 
the previous study) . The results of these simulations are 
resumed in Fig.|3 In this figure we also show the output 
of some 2D simulations, which confirm the observation 
that ID simulations provide a good estimate for the lo- 
cation of the spiral-antispiral transition. Also here it is 
seen that the prediction based in the CGLE description 
is valid only in the vicinity of the Hopf threshold. As 
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FIG. 5: Plot of the measured wavenumber (js selected by ID- 
and 2D-(anti)spirals (open circles and squares, respectively) 
in numerical simulations of the FHN model (IIH as a function 
of the parameter e, with \/5 = 1 and /3 = (in this case 
7 = 1/2 and e'n = 2). The full lines indicate the wavenumber 
selected by the 2D antispirals in and ID antispiral analogues 
in the CGLE. The CGLE ID line is calculated using the 
wavenumber q^i predicted analytically by Hagan's formula 
0. In the inset the value of g| in the region close to the 
onset is shown. Note that the CGLE ID and CGLE 2D 
lines are directly proportional to the distance to the threshold 
^ejj- — e and consequently that both qsi and qs are constant 
in the framework of the amplitude equations, but also note 
that this constant value differs [3^ . The measured g^i and 
qs approach (as expected) the predictions of the ID and 2D 
CGLE antispirals near the threshold but are systematically 
smaller far away from it. 



the distance to this threshold increases the antispirals 
arc turned into spirals. 



IV. SUMMARY 

In this paper we emphasise the fact that the essen- 
tial difference between spirals and antispirals is the sign 
of phase velocity of the travelling waves emanating from 
the core region and that phase and group velocity do not 
necessarily have to point into the same direction. Con- 
sequently, we may conclude that spirals and antispirals 
can not occur simultaneously. They appear in different 
regions of the parameter-space. 

We have shown that spirals and antispirals are related 
by a simple symmetry transformation in the parameter 
space of the CGLE and derived a criterion for the CGLE 
coefficients, namely ci + > 0, that predicts the ex- 
istence of antispirals in the respective RD model near 
a Hopf instability. In the CGLE parameter space, an- 
tispirals appears as frequently as spirals. Spirals and 




FIG. 6: Summary of the results of numerical simulations 
of ID- and 2D-(anti)spirals in the Brusselator model with 
\/5 = 1.5. The thick full line indicates the location of the Hopf 
instability threshold and the white circles indicate the spiral- 
antispiral transition points at onset (c/. Fig. El . The dashed 
curves indicate (approximately) the analogue of the spiral- 
antispiral transition in the ID calculations. The squares and 
diamonds show the location of some 2D numerical simulations 
resulting in spirals and antispirals respectively. 



antispirals arc found in complementary regions of the 
paramtcr space and arc characterized by the selection of 
positive and negative wavenumbers, respectively. At the 
boundary between antispirals and spirals, the selected 
wavenumber is zero and the selected wavelength diverges. 
Our criterion has been mapped to the parameter space of 
two representative reaction-diffusion models. The anal- 
ysis of both models suggests that to get antispirals it is 
necessary that the inhibitor also diffuses. Moreover, the 
diffusion constants of both reactants should be of simi- 
lar magnitude in order to find antispirals, which is likely 
to occur near the codimension-2 Turing-Hopf instabil- 
ity. We have also performed numerical simulations of 
both models near and far form the oscillatory instability 
threshold. Antispirals arc quite likely to appear near the 
Hopf bifurcation but turn into regular spirals far into the 
oscillatory where the amplitude of waves becomes large. 
Both examples show that the assumption that the ID 
spiral analog gives a good approximation of the spiral- 
antispiral transition behaviour not only near onset but 
also moderately far from it. 

In a more general distinction, antispirals emit phase 
waves that are typical for oscillatory media. In contrast, 
trigger waves common in excitable media should always 
have the same sign of group and phase velocity. Since 
the distinction between trigger and phase waves is not 
sharp, we cannot make a general statement about the 
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possibility of antispirals in excitable media, but expect 
that they will be quite hard to find therein. The rarity 
of experimental observations of antispirals suggests that 
most experiments exhibiting spirals are conducted either 
far away from the oscillatory threshold or under excitable 
resp. bistable conditions. It is tempting to explain the 
experiment of Vanag and Epstein [Tl| from the proximity 
of the control parameters to a Hopf bifurcation, though 
other more complicated possibilities cannot be ruled out 
at this stage. 

Finally, we like to stress the most essential point in 
the consideration of the nature of spirals and antispirals. 
The crucial quantity for both structures is the group ve- 
locity of the periodic waves outside the core region. In all 
models considered, the group velocity is pointing outward 
from the spiral core as has been shown analytically for 
the CGLE. The distinction between spirals and antispi- 
rals stems only from the sign of the phase velocities and 
the corresponding phcnomcnological impression. In our 
opinion is it important not to confuse the distinction be- 
tween spirals and antispirals with the distinction between 
sinks and sources. Several authors have pointed out that 
the definition of a source depends only on the group ve- 
locities and not at all on the phase velocity of the periodic 
waves far away from it [40, HH . Counting arguments show 
that sources belong typically to a zero-parameter family 
while sinks are members of a two parameter family [40j . 
In other words, sources select frequency and wavelength 
of the emitted periodic waves, while sinks simply appear 
where two wavetrains of arbitrary wavelength and period 
happen to meet. A classification of defects in reaction- 
diffusion media can be done according to the group veloc- 



ities of periodic waves on both sides of the defect and the 
velocity of the defect itself. Apart from sources and sinks, 
one may also find transmission defects (one-parameter 
family) and contact defects in reaction-diffusion media 
|4ll |. In the present context, spirals and antispirals are 
both sources, i. e. they select their frequency and wave- 
length and control their surroundings. If they were sinks 
(as sometimes claimed in the literature) , the waves had to 
be organized from the domain boundaries, which would 
require another unknown "organizer" structure. Alto- 
gether, we have shown by combination of existing ana- 
lytical arguments that antispirals are structures of sim- 
ilar nature as spirals and should hence be expected in 
many systems exhibiting phase waves in particular near 
the onset of oscillations. Their experimental discovery 
can be nicely interpreted within the existing theoretical 
framework for oscillatory reaction-diffusion media. 
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